The effect of vacancy-induced magnetism on electronic transport in armchair carbon nanotubes 
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The influence of local magnetic moment formation around three kinds of vacancies on the electron con- 
duction through metallic single-wall carbon nanotubes is studied by use of the Landauer formalism within the 
coherent regime. The method is based on the single-band tight-binding Hamiltonian, a surface Green's func- 
tion calculation, and the mean-field Hubbard model. The numerical results show that the electronic transport 
is spin-polarized due to the localized magnetic moments and it is strongly dependent on the geometry of the 
vacancies. For all kinds of vacancies, by including the effects of local magnetic moments, the electron scat- 
tering increases with respect to the nonmagnetic vacancies case and hence, the current-voltage characteristic of 
the system changes. In addition, a high value for the electron-spin polarization can be obtained by applying a 
suitable gate voltage. 

PACS numbers: 72.10.-d, 72.10.Fk 



INTRODUCTION 



transport especially near the Fermi energy. 



The electrical transport properties of single-wall carbon 
nanotubes (CNTs) and other carbon-based materials have at- 
tracted much attention due to their unusual properties and 
great potential for technological applications |ll|-|4j]. Among 
these features, the ballistic electron conduction and the long 
range spin coherent transport for perfect and defective single- 
wall CNTs have been investigated theoretically and experi- 
mentally iil-lTtl- These important properties of spin polarized 
electrons in single-wall CNTs has motivated their use in the 
emerging field of spin electronics |8|] which aims to effec- 
tively control and manipulate the spin degrees of freedom in 
the electronic devices 

In order to achieve spintronic devices based on single-wall 
CNTs, it is important to understand the magnetic effect of va- 
cancies and impurities on the electron conduction (js, 11|. 
because, the electronic properties of carbon-based materials 
strongly depend on their topological structure 11211 . Therefore 
the electronic structure of the CNTs can differ due to the topo- 
logical defects or the addition of different compounds. On 
the other hand, in the transport processes, the ballistic con- 
ductance depends on the number of conducting channels at 
the Fermi energy 11311 . Consequently, the appearance of the 
vacancy defects in a structure can change the electronic and 
transport properties of the system jstl- Also, the conductance 
of an imperfect system is lowered due to the reflection of elec- 
tron waves from the defects iflill . In addition, the localized 
states near the vacancy are magnetic and change the net mag- 
netic moments in the carbon-based nanostructures. The hon- 
eycomb lattice of graphene sheet is formed by two sublattices 
A and B (bipartite lattice). For a bipartite lattice with differ- 
ent numbers Na and Nb of sites and the Hubbard repulsive 
parameter, the total spin, S, of the ground state of the system, 
which is mainly localized near the vacancy, is 25* = Na — Nb 
1 1511 . This important feature can block and change the spin 



In this regard, the electronic and magnetic properties of 
vacancies m smg le-wall CNTs d and graphene nanorib- 
bons II17II have been investigated. In addition, the electron- 
spin polarization has been observed in the CNTs when doped 
with magnetic adatoms or molecules lIlSll . Moreover, the lo- 
calized states of the impurity can change the spin-polarized 
conduction in the presence of a gate voltage or applied bias 
lUKl^. Recently, the effect of vacancy on the conductance of 
smg le-wafl CNTs M, 

I2III and the spin-dependent transport 
properties in ferromagnetically contacted single-wall CNTs 
have been investigated | 22 , 23tlj but the magnetic behavior 
of vacancy defect as regards the transmission of single-wall 
CNTs has not been considered. 



The purpose of this work is to study the effect of magnetic 
vacancies on the spin-polarized transport through armchair 
CNT junctions and manipulate this polarization by means of 
gate voltages. We simulate ideal vacancies which are made 
by removing carbon atoms from lattice sites without includ- 
ing the deformation of the tube wall around the vacancies. We 
consider three typical vacancy types. For the first type, a sin- 
gle carbon site (A or B) is removed. For the second one, two 
same sites (two A or B sites), and in the last type, four A sites 
in one carbon ring are removed, as shown in Figs. 1(a) and 
1(b), respectively. 

Using the single-band tight-binding approximation and the 
mean-field Hubbard model ll24ll . the electronic structure and 
the localized magnetic moments around the above-mentioned 
vacancies are calculated. Also, using the non-equilibrium 
Green's function technique and the Landauer-Biittiker the- 
ory 11251], the spin-polarized transport in defective single-wall 
CNTs is investigated. 
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FIG. 1: Five unit cells of a (4,4) single-wall CNT (consist of 80 
atoms) as the channel, in the presence of (a) two and (b) four vacan- 
cies on a ring of carbon atoms. 



MODEL AND FORMALISM 

We consider a system that consists of a defective CNT with 
finite length, sandwiched between two perfect semi-infinite 
CNTs as electrodes. All CNTs have same structure, (n,n). To 
calculate the electron transmission properties of this defective 
CNT which acts as a channel, the electronic structure of this 
system should be resolved in detail. Hence, we decompose 
the total Hamiltonian of the system as ll26ll 



(1) 



where Hl and Hj^ are the Hamiltonians of the left (L) and 
right (i?) electrodes. He describes the channel Hamiltonian 
and contains the magnetic properties of the vacancies, and 
Ht is the coupling between the electrodes and the central part 
(defective CNT). The total Hamiltonian is described within 
the single-band tight-binding approximation. Therefore, the 
Hamiltonian of the electrodes can be written as ll27ll 



t- - )c 



(2) 



<«cjcv>,0- 



where c| ^ (ci^^a) creates (destroys) an electron at site i in 
electrode a(= R,L) and the hopping elements ti^j^ are 
equal to t for nearest neighbors and zero otherwise. Here, ea 
is the on-site energy of the electrodes and will be set to zero. 
The coupling Hamiltonian is described as 

Ht = - E *'-^o(cL,^^^c,-+H.c.) , (3) 



The parameters for hopping between the electrodes 

and the channel are taken to be t. In order to obtain the mag- 
netic moment of each atom around a vacancy, we use the Hub- 
bard model within the unrestricted Hartree-Fock approxima- 
tion 1241]. Accordingly, the Hamiltonian of the central part 
within the mean-field approximation of the Hubbard model 
can be written as 12412811 



+ {nii)nti - {nit){nii)} , (4) 

where is the on-site energy and will be set to the gate po- 
tential in the channel region, the second term corresponds 
to the single 7r-orbital tight-binding Hamiltonian, while the 
third term accounts for the on-site Coulomb interaction U. In 
this expression, dl^{dia-) creates (annihilates) an electron, and 
{fiia-) is the expectation value of the number operator for an 
electron with spin a at the ith site. 

The Green's function of the system can be written as 



Gcie) = [ei -He- Sl(£) - ±Ris)Y 



(5) 



where ^^(e) = Hca9a{£)Hlj^ is the self-energy matrix 
which contains the influence of the electronic structure of 
the semi-infinite electrodes through the lead's surface Green's 
function Qa lEgilsoll . Correspondingly, He a defines the ma- 
trix of coupling between the surface atomic orbitals of the 
lead a and the channel. We should note that, Gc is a spin- 
dependent matrix of size 2N, with N being the number of 
atoms of the central region. Therefore, due to our mean-field 
decoupling scheme, we can decouple the electronic transport 
into spin-up and spin-down currents. 

Using the iterative method for the calculation of the trans- 
fer matrices based on the principal layers concept, the surface 
Green's function of the left and right leads can be obtained as 

lii: 



gL{e) = [ei-H^,,-H^^fL]. 
gnie) = [ei - H^o - H^iTr] , 



(6) 



where Hqq and Hq^ are matrices and correspond to an iso- 
lated principal layer and the interaction between two nearest 
principle layers in lead a, respectively. Also, Tl and Tf> are 
the transfer matrices and for the details of the calculations the 
reader can refer to Ref ll29ll . From the expression of the non- 
equilibrium Green's function, the spin-dependent local den- 
sity of states (LDOS) and the expectation value for the number 
operator of electrons in the channel are given by 



D,. 



1 



lm{{ia\Gc{e)\ia)) 



(7) 
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Accordingly, the magnetic moment, m^, at site i of the chan- 
nel can be calculated using 



it) 



(8) 



In this study, we solve the mean field Hamiltonian self- 
consistently by an iterative method 12411 . In the first step, 
we start from an anti-ferromagnetic configuration as an ini- 
tial condition and establish the Hamiltonian for the Hubbard 
model [Eq. ^]. In the second step, the effect of electrodes 
on the channel is added via and E/j, and then the Green's 
function of the channel, Gc, in the presence of the electrodes 
is calculated. In the third step, the expectation values of the 
number operators (riicr) at each site and for both spin align- 
ments are calculated by using the Green's function. Finally 
the new expectation values of the number operators are re- 
placed in Eq. (01), and this process is repeated until the dif- 
ference between two successive iterations becomes less than 
10-1 

The total Hamiltonian does not contain inelastic scattering, 
because there is no spin-flip or the other sources of scattering 
in the system. In other words, the transmission probabilities of 
majority and minority subbands are independent and the elec- 
tronic transport can be decoupled into two spin currents: one 
for spin-up and the other for spin-down. Therefore, the spin- 
dependent currents for a constant bias voltage are calculated 



by using the Landauer-Biittiker formalism Ill0ll26ll 



e 



)^ fie- ^lR)]de , (9) 



where / is the Fermi-Dirac distribution function, hl.r = 
Ep ± \eVa are the chemical potentials of the electrodes 
and Tcr(e) = "^t:^ lG c r&q <j\ energy- and spin- 

dependent transmission function. Using S^, the coupling ma- 
trices Ta can be expressed as Fq = — 2Im[SQ(£)]. 

RESULTS AND DISCUSSION 

We study the LDOS, the transmission and the current- 
voltage characteristic of the single-wall CNTs in the presence 
of vacancies when the CNT in all regions is of (4,4) type. We 
expect the magnetic behavior of the vacancy to affect the elec- 
tron conduction through the structure, since different scatter- 
ings occur for the electronic waves with different spin densi- 
ties. Therefore, the transmission for majority electrons can be 
different from that of the minority electrons at certain ener- 
gies. 

In the cases of the first, second and third vacancy types, the 
total magnetic moments for the finite CNT (channel) as shown 
in Fig. 1, are 1, 2 and 4/_tB respectively, which are in good 
agreement with the Lieb's theorem llisll . However, when the 
two semi-infinite CNTs (electrodes) are attached to the chan- 
nel, their electronic properties are changed. This causes ap- 
preciable variation in the value of magnetic moment of each 
atomic site near the vacancies and affects the spin-dependent 




FIG. 2: LDOS per site around two vacancies for the CNT/defective- 
CNT/CNT junction with U = 1.06 1, Vg = 0.0 and Va = 0.0 1. 
Full (open) circles are for majority (minority) spin electron. The dot- 
ted (dashed-dotted) line is for majority (minority) spin in the absence 
of CNT leads ( and the dashed line indicates the Fermi energy). Note 
that, the majority and minority electrons are determined by Eq. Q. 




FIG. 3: Transmission coefficients as a function of energy in the pres- 
ence of two vacancies (second type) with Va = l.Qt and Va = 0.0 1. 
Full (open) circles are for majority (minority) spin electrons. 



transport through the channel. When the atoms from the same 
sublattice are removed, the final spin configuration, after self- 
consistent calculations for the atomic sites near the vacancies, 
is ferromagnetic. These magnetic moments only obtain for 
the atomic sites near the vacancies and for the other sites ap- 
proximately vanish. Our numerical results showed that in the 
presence of a single vacancy, the magnetic moment vanishes 
approximately and an unpolarized ground state is obtained for 
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the value of J7=1.06 t, due to the effect of semi-infinite elec- 
trodes on the channel 131 1. 

In order to clarify this effect, we have plotted in Fig. 2 
the LDOS in the presence of two vacancies. In a defective 
CNT, the electronic states on the carbon atoms near the va- 
cancy strongly depend on the presence or absence of CNT 
electrodes. In a finite-length CNT, the electronic states on 
carbon atoms are almost localized due to the absence of pe- 
riodicity along the axis of CNT. Therefore, sharp features are 
observed in the spin-dependent LDOS for the states around 
the Fermi energy. When the electrodes are attached to the de- 
fective CNT, the localization of electronic states is suppressed 
and the sharp peaks are broadened for both majority and mi- 
nority electrons. Accordingly, the difference between the two 
densities of states and hence the spin polarizations are re- 
duced. Note that in both cases the condition (fii-^) + {ui^^ ~ 1 
at each site is satisfied. 

It should be pointed out that the Hubbard repulsive param- 
eter shifts the Fermi energy from zero to Ep = 0.3 1 as shown 
with a vertical line in Fig. 2. In fact, when the two vacan- 
cies in the channel are included in the calculations, we ob- 
tain 5'=0.55 for the total spin value in the presence of CNT 
electrodes and the maximum value of the magnetic moment 
reaches 0.16 /is. This magnetic moment can affect the spin 
transport in the channel. 

To see the effects of vacancy-induced magnetism on the 
electron conduction, we have shown in Fig. 3 the transmis- 
sion coefficients as a function of energy for majority and mi- 
nority electrons at VG=l-Ot. The transmission spectrum of 
two spin subbands is non-degenerate for all electronic states 
within the energy window and there is a significant difference 
between the two spin subbands. Since in the process of elec- 
tronic transport, the electrons with Fermi energy carry most 
of the current, we have shown the transmission coefficients 
around the Fermi level. It is clear that the Fermi energy shifts 
from 0.3 1 to 1.3 1 due to the gate voltage. 

We note that a gate voltage shifts the electronic states and 
hence, the transmission for two spin subbands in the central 
region can vary significantly. Therefore, a gate potential is 
able to change the spin transport through such a defective 
CNT. In order to see the effect of the gate voltage on the spin 
transport, we have calculated the degree of spin polarization 
(P) for electrons traversing the channel which can be defined 



as P{E) 



Fig. 4 shows this physical quantity 



in terms of the gate voltage in the case of two vacancies. It 
is clear that the spin polarization can reach values as high as 
65%. Also, we found that the change in the position of the 
two vacancies with respect to each other causes appreciable 
variations in the value of the localized magnetic moment and 
hence, in the degree of spin-polarization. Generally, on in- 
creasing the distance between two vacancies the total value of 
the channel magnetic moment decreases. 

Fig. 5. shows the transmission coefficients as a function 
of energy for majority and minority electrons in the case of 
four vacancies (third type). Only the energy window around 
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FIG. 4: Degree of spin polarization as a function of gate voltage for 
vacancies of the second type at the Fermi energy and Va ~ 0.0 1. 
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FIG. 5: Transmission coefficients as a function of energy in the pres- 
ence of four vacancies (third type), for Vg ~ 1-0 1 and Va = 0.0 1. 
Full (open) circles are for majority (minority) spin electrons. 



the Fermi level has been shown {Ep ~ 1.3 1). As we know, 
in the present structure where the electrodes are made of 
CNTs, there are more paths for the electrons to pass from the 
electrodes to the defective CNT. Therefore, in some energy 
ranges of this figure, the spin-dependent transmission coeffi- 
cients can be larger than unity. The total value of the spin is 
1.40 and the maximum value of the magnetic moment reaches 
0.30 HE- Also the maximum value of the spin polarization at 
Vq = 1.0 1, reaches 80% which indicates that the local mag- 
netism induced by vacancies can be useful in spintronic de- 
vices. This spin polarization of the electrical conductivity can 
play an important role in the current- voltage characteristics in 
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the channel. 

For this reason and to further understand the effect of 
vacancy-induced magnetism, we have plotted in Fig. 6 the 
spin currents in terms of the applied voltage. Since the elec- 
tron conduction depends on the electronic states lying be- 
tween fiL and flu, on increasing the applied voltage, these 
states move inside the energy window and the electronic cur- 
rents for both majority and minority spin electrons increase. 
From the current-voltage characteristic, we see that the ap- 
plied voltage can change the difference between both spin-up 
and spin-down currents significantly and therefore, we should 
adjust the value of the gate and applied voltages to obtain max- 
imum values for the degree of polarization of the spin cur- 
rents. 

In this study, we only focused on the A-site defects because 
according to Lieb's theorem, when the same sublattice atoms 
are removed the total magnetic moment increases and there- 
fore the polarization in the spin currents becomes remarkable. 
Also, our numerical results showed that if we remove two dif- 
ferent sites from one carbon ring, the local magnetic moments 
near the vacancies are nonzero, while the total magnetic mo- 
ment is zero and hence, the transmission coefficients of two 
spin channels are degenerate and practically unpolarized. In 
other cases, i.e., when different numbers of A and B sites are 
removed, the total magnetic moment is small and a little spin 
polarization can be observed (not shown here). 
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FIG. 6: Cuffent- voltage characteristics in the presence of four vacan- 
cies. The solid (dashed) line is for majority (minority) spin electrons. 
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CONCLUSION 

In this work, the coherent spin-polarized transport through 
a CNT/defective-CNT/CNT junction was investigated on the 
basis of the non-equilibrium Green's function technique, and 
the mean-field Hubbard model. We have shown that there are 
major differences in the behavior of the electronic transport 
when different kinds of vacancies are considered. The nu- 
merical results indicate that, the point defects of which dif- 
ferent types were considered-and hence with different mag- 
netic responses-could effectively alter the electron conduction 
through the junction. 

In the defective CNTs, the difference between the transmis- 
sion coefficients of two spin subbands at some energy values 
is small. Consequently, in such a case, the degree of spin po- 
larization can be manipulated by applying a gate voltage to the 
channel. In addition, our results suggest that by choosing the 
position of the vacancy in the channel properly to enhance the 
local magnetization, one can increase the values of the spin- 
polarized currents. These defects may be useful in the process 
of spin injection into semiconducting devices for spintronic 
applications. 

Throughout the study, we have ignored the effect of spin- 
flip scattering, especially near the vacancies. This factor may 
affect the spin-dependent transport, and hence, a further im- 
proved approach is needed to obtain more accurate results. 
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